Violation of Onsager Reciprocity in Underdoped Cuprates ? 
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One of the canons of condensed matter physics is the Onsager Reciprocity principle p] in systems 
in which the Hamiltonian commutes with the time-reversal operator. Recent results of measurements 
of the Nernst coefficient in underdoped YBa2Cu30e+ x , together with the measurements of the 
anisotropy of conductivity and the inferred anisotropy of the thermopower, imply that this principle 
is violated. The probable violation and its temperature dependence are shown to be consistent with 
the Loop-current phase which has been directly observed in other experiments. The violation is 
related directly to the magneto-electric symmetry of such a phase in which an applied electric field 
generates an effective magnetic field at right angle to it and to the order parameter vector, and 
vice-versa. 



Extraordinary Experimental results 2J for magneto- 
transport experiments have been recently presented for 
various dopings in YBCO. Starting below the pseudogap 
temperature T*(x), the magnitude of the Nernst coeffi- 
cient rises. Persuasive evidence has been presented that 
except for a narrow region ST C /T C < 1/4, it is due to 
quasi-particles. We are not concerned here with its over- 
all magnitude, but only with its symmetry. Nernst effect 
is the chiral response of fermions to mutually orthog- 
onal magnetic field and a thermal gradient. We focus 
attention on the observation that the observed chirality 
in Nernst effect is non-reciprocal (NRN) : For a magnetic 
field in the z— direction, the absolute magnitude of its 
value \v a b\ is found [2] to be different with the thermal 
gradient applied in the b direction and voltage measured 
in the a direction, from \v ba \, the value with the thermal 
gradient applied in the a direction and voltage measured 
in the b direction, a and b are the directions of the or- 
thorhombic crystalline axes in the Cu-0 planes, b is taken 
to be along the chains. A suitable definition of the NRN 
in this case, which avoids issues of sign of the Hall effect, 
thermopower etc., which are not relevant to the present 
discussion, is 



NRN : 
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\Vab\ 
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In Ref. ©, it is concluded that NRN varies smoothly 
above T*(x) but that it has a contribution below T*(x) 
consistent with ex (T*(x) - T) for small (T*(x) - 
T)/T*(x), saturating to a constant value at lower tem- 
peratures, i.e. it has a contribution similar to the square 
of an order parameter in the mean- field regime of a phase 
with transition temperature T*{x). This behavior is ob- 
served in the range 0.1 < x < 0.18. In these experiments, 
as well as in others, the uncertainty in determination of 
T*(x) is ps ±20K. We may isolate the temperature de- 
pendent contribution below T* (x) by defining 

X{T) = NRN(T) - NRN(T ^ T*(xj), T < T*(ar).(2) 

Related quantities are deduced in Ref. ©. The Nernst 



coefficient, linear in B, under the condition of zero- 
transverse electrical current is given by [5] 



Vba 



O aa&ba &ba&aa 
Vaa&bb 



(3) 



In Eq. J3j a and a are the conductivity and thermoelec- 
tric tensors. For v ab , the subscripts, (a, b) —> (b,a). 
There are two classes of possibilities for NRN: 

(I) Diagonal anisotropy only : a aa ^ a bb or a aa ^ a bb or 
both. 

(II) Off-diagonal anisotropies a ab =^ — a ba or a ab ^ —o ba 
or both. This may or may not be accompanied by diag- 
onal anisotropies. 

YBCO is orthorhombic, so anisotropies in the diagonal 
components of a and a are to be expected at all temper- 
atures. The anisotropy has been explicitly measured for 
the conductivity [2HI]; it is weakly temperature depen- 
dent and varies quite smoothly across T*(x). Explicit 
results for the anisotropy of thermopower are not avail- 
able for many compositions. The ratio measured for the 
Oe.35 compound is slowly varying at all temperatures [5]. 
Thermopower measured [1] along the a-direction at vari- 
ous dopings of interest is just as smooth across T*(x) as 
a aa and quite unlike the Nernst coefficient. In general, it 
is improbable that thermopower show non-analytic be- 
havior while the conductivity does not. This view is 
strengthened by the available data [H [5] . Therefore, pos- 
sibility / can be excluded for X(T). 

Consider possibility II. If the states of the system are 
eigenstates of the time-reversal operator Onsager Reci- 
procity relations [U |S] hold so that a ab = —a ba for any 
crystalline symmetry. Also, if as is usual the energy re- 
laxation rate is related to the momentum relaxation rate 
so that the Mott relation is obeyed, a ab — —a ba . So for 
a state with time-reversal preserved NR must be smooth 
across T ps T* (x) . The experimental observations there- 
fore imply possible time-reversal violation in the pseudo- 
gap phase. However as seen below, this is only necessary, 
not sufficient; not all time-reversal violating states and 
crystalline symmetries lead to the observed behavior of 
X(T). For example, ferromagnetism or states with sim- 




FIG. 1: One of the domains of the observed order. A rect- 
angular unit-cell is shown with orbital-magnetic moments in 
the + and - z directions in two of the four O-Cu-0 triangles. 
The anapole order parameter L for this domain is also shown. 
Note that according to p}, L in general does not point along 
the diagonal of the rectangle. 



pie spin or orbital antiferromagnetic symmetries cannot 
give the observed behavior. The experiments may then 
be a macroscopic manifestation of time-reversal breaking 
with rather special features. 

Transition to states with time-reversal violation (TRV) 
were predicted with an onset temperature T* (x) [7J with 
symmetries consistent with observations in four differ- 
ent families of cuprates, by polarized neutron scattering 
[8] or dichroic ARPES 9. Polar Kerr effect consistent 
with fcrromagnetism with a small moment has also been 
observed [TO] and is discussed in Ref. [TT]. We show 
that such a state in an orthorhombic crystalline symme- 
try does indeed give a NR in the Nernst effect consistent 
with the experiments. 

The symmetry of the observed TRV state has been dis- 
cussed in Ref. [Jj [12] . The order parameter violates time- 
reversal and reflection symmetries but preserves transla- 
tional symmetry and is given by the polar time-reversal 
odd vector, an anapole- vector, 



d 2 r (m(r) x f). 



(4) 



unit — cell 



m's are the generated orbital magnetizations due to 
electron-electron interactions below T*(x). They point 
along ±z in two of the opposite O-Cu-0 triangles in 
the unit-cell. One of the four possible domains in a or- 
thorhombic crystal with m and L displayed is shown in 
Fig. (fl]). In an orthorhombic crystal, the four domains 
have 

L = ±L x x±L y y. (5) 

We will take x, y to be along the a and b— directions of 



the crystal respectively. In a tetragonal crystal L x = L v , 
but not in an orthorhombic crystal. 

As has been noted [7J H2] this state has Magneto- 
electric symmetry |13| : An external field Bz, generates 
an electric field, or more generally any vector r of the 
symmetry of an electric field (polar- time-reversal even) 
in the plane; an external electric field E in the plane 
generates a magnetic field, or more generally any vector 
(r x v) of the same symmetry as a magnetic field (axial- 
time-reversal odd) along z. Such generated vectors in a 
metal can be represented by fermion operators, so that 
for L ^ 0, the Hamiltonian acquires terms, 



H' = x((Lx B)^ f (r)er^(r) 

, j. . ^ e(r x v — v x r) 



(6) 



-(LxE).f(r) 



2c* 
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if), ip + are the fermion annihilation and creation operators 
and x is the magneto-electric coefficient. We have cho- 
sen units such that electric and magnetic fields have the 
same dimensions so that \L is dimensionless. We have 
also introduced a phenomenological velocity c* , which or- 
dinarily would be the velocity of light c. However, since 
the effective electric and magnetic fields are generated 
through electron-electron interactions, it is quite likely 
that c* is related to the microscopic parameters of the 
model and so much less than c. Only a microscopic cal- 
culation of the magneto-electric coefficient can tell; at 
the moment we leave it as a parameter to be determined. 
The force on the fermions due to H' is obtained by cal- 
culating the commutator with the momentum operator 
p to get 



dp ( , „ N (L x E) x v 



(7) 



We see that in addition to the external electric E and 
magnetic B fields, internal electric e and magnetic b 
fields act on the fermions: 



e = x (LxB), b=- X (c/c*)(LxE). 



(8) 



Standard transport theory [14] now gives that charge- 
Hall current j^ is given by 



;H 



cr H (B + b) x (E + e) 



(9) 



where o~ H is the normal Hall conductivity tensor. For- 
mally, reciprocity is obeyed only for interchange of direc- 
tion of charge-Hall current and of (E + e). Looking at 
([8]), this cannot in general be achieved in experiments. 
Similarly the thermal Hall current is given by 



Jt 



a H (B + b) x (E + e) 



(10) 



where a H is the normal thermal Hall conductivity tensor 
which in the present case formally obeys reciprocity only 
for interchange of direction of thermal-Hall current and 
of (E + e). 



Eqs. ([9]) and ( 10 ) have the usual terms proportional to 
E x H, which do not contribute to the NR defined in Eq. 
(fTl). From the cross-terms in E and b and in B and e, 
one gets respectively Hall effects with no external B or 
no external E; they are linear in L and therefore cancel 
in transport measurements in macroscopic samples in the 
sum over different domains of L. The terms of immediate 
interest are proportional to 



exb = x 2 (c/c*)(L-(BxE))L. 



(11) 



Unlike the conventional Hall current, this contribution 
is " uniaxial" ; it flows only along L and depends on the 
projection of B x E onto L. This cannot obey Onsager 
reciprocity. The maximum NR is for the electric field ap- 
plied parallel or perpendicular to L, where the anomalous 
contribution is zero for the latter case. 
On using Q and Q give contributions, 



o ab = a H X 2 (c/c*)Ll 
a ba = a H X 2 (c/c*)Ll 



(12) 



These contributions are the same for all four domains. 
There are similar non-reciprocal contributions in a ab and 
ata- 

To calculate the value of NR, we need both the off- 
diagonal and the diagonal components of the transport 
tensors. <7 a a/cbb is measured [21 0]. As discussed above, 
based on available data we assume that a aa /abb behaves 
similarly. Thus, 



0~bb 



abb 



(13) 



with a nearly temperature independent r in the range of 
interest. The measured value OSj is about 1/2. 

Then for T just below T*(x) where (c/c*)x 2 L 2 << 1 



X{T)^(c/c*) X \Ll-Ll) 



(14) 



This agrees with the experimental result that X(T) « 
(T*(x) — T), in the mean-field range in which L 2 oc 



(T*(x)—T). At low temperatures, where L 2 is a constant, 
X(T) is also a constant. This is also consistent with the 
experiments. The sign of the effect is also consistent if 
L 2 > L 2 , as is expected. 

Given the available data, it is impossible to calculate 
the magnitude of the effect and compare with the exper- 
iments. The problem is (1) that since r is neither too 
small compared to 1 or too close to 1, the expressions for 
X(T) cannot be simply disentangled in terms of L and 
r without separately knowing the magnitudes of a a b and 
a a b for at least some temperature, and (2) r has a con- 
tribution due to chains which cannot be disentangled in 
the experiment so that a contribution to the magnitude of 
the Nernst coefficients from the planes can be separately 
obtained. It would be far more satisfactory to directly 
measure the ratio of the off-diagonal parts, o~ a b/o~ba and 
a a b/ab a to test for violation of Onsager reciprocity di- 
rectly, unlike the Nernst coefficients which has diagonal 
and off-diagonal parts of the two different transport ten- 
sors mixed in. We urge such experiments. Especially 
interesting would be the angular dependence suggested 
through Eq. (fTT|. 



To summarize, the experimental results on the 
anisotropy of the Nernst coefficient are consistent with 
a violation of the Onsager Reciprocity for T < T* (x) if, 
as reported, the diagonal components of the conductiv- 
ity, and as is probable, the diagonal components of the 
thermopower are smooth across T*(x). Such a violation 
is consistent with the magneto-electric symmetry already 
deduced in microscopic experiments in Cuprates. 
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